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Abstract: We discuss connections between certain classes of supersymmetric quiver gauge 
theories and integrable lattice models from the point of view of topological quantum field 
theories (TQFTs). The relevant classes include 4d AA = 1 theories known as brane box 
and brane tilling models, 3d M = 2 and 2d J\f = (2, 2) theories obtained from them 
by compactification, and 2d A/" = (0, 2) theories closely related to these theories. We 
argue that their supersymmetric indices carry structures of TQFTs equipped with line 
operators, and as a consequence, are equal to the partition functions of lattice models. 
The integrability of these models follows from the existence of extra dimension in the 
TQFTs, which emerges after the theories are embedded in M-theory. The Yang-Baxter 
equation expresses the invariance of supersymmetric indices under Seiberg duality and its 
lower-dimensional analogs. 
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1 Introduction 

The present work is motivated by an intriguing connection discovered in the past few years 
between supersymmetric quiver gauge theories and integrable lattice models in statisti¬ 
cal mechanics. 

In [1, 2], Bazhanov and Sergeev introduced an integrable spin model on a planar 
lattice, which generalizes many of previously known integrable lattice models such as the 
Kashiwara-Miwa [3] and chiral Potts [4-6] models. Soon after the first paper by Bazhanov 
and Sergeev appeared, Spiridonov [7] gave an interpretation of their model in terms of 
Af = 1 quiver gauge theories in four dimensions. The relation to gauge theory was further 
elucidated by Yamazaki [8, 9], who realized that the relevant quiver gauge theories arise 
naturally from a particular class of brane configurations in string theory, called brane 
tilings [10]. From the gauge theory viewpoint, the lattice on which the spin model is 
defined is the quiver diagram, the partition function is the supersymmetric index, and the 
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Yang-Baxter equation that guarantees the integrability of the model is a special instance 
of Seiberg duality. 

This discovery, while quite remarkable, leaves us with a series of questions: Why are the 
supersymmetric indices of these theories captured by a lattice model? Why is this model 
integrable? Are there structures of integrable lattice models hidden in other theories? 

In this paper we answer these questions, combining ideas from two equally stimulating 
connections uncovered in recent years. Of these connections, one lies between a certain 
class of quiver gauge theories and topological quantum field theories (TQFTs), whereas 
the other relates TQFTs and integrable lattice models. Our goal is to connect quiver 
gauge theories and integrable lattice models, and TQFTs provide a bridge between the 
two elements. 

The first connection in question arises from the M5-brane construction of 4d A? = 2 
theories [11-13]. Consider a stack of M5-branes wrapped on x x S, where S is a 
compact Riemann surface. In addition, we introduce M5-branes that intersect with these 
branes along submanifolds of the form x x {pi}, with pi being points on S. If S 
is small compared to the and S^, this brane system is described at low energies by 
an A? = 2 theory on x S^. Often this is a quiver gauge theory. The path integral 
on the geometry x (after analytic continuation to Euclidean spacetime) computes 
the supersymmetric index of the theory. The most important property of this quantity 
is that it is invariant under continuous changes of the parameters of the theory. Since 
the geometry of S encodes such parameters, it follows that the supersymmetric index is 
a topological invariant of S. In fact, it is equal to a correlation function in a 

TQFT on S, where Oi is a local operator representing the M5-brane inserted at pi [14]. 
If we consider a protected quantity different from the index on (such as the lens space 
index [15]), then we get a different TQFT. 

The second connection refers to a general construction of integrable lattice models 
from TQFTs of a special kind, due to Costello [16, 17]. Given a 2d TQFT equipped 
with line operators, we can place it on a torus and wrap line operators around various 
1-cycles that form a lattice. The topological invariance of the theory implies that the 
correlation function for this lattice of line operators coincides with the partition function 
of a statistical mechanics model defined on the same lattice. Under this correspondence, 
the Yang-Baxter equation translates to the statement that the correlator remains the same 
when a line operator is moved past the intersection of two other line operators. The 
structure of a TQFT itself is not strong enough to ensure this property since the move 
is not topologically trivial; the system may undergo a phase transition as the moved line 
hits and crosses the intersection point. If, however, the theory has “extra dimensions” 
along which the line can be shifted, the collision can be avoided and hence the Yang- 
Baxter equation holds. What is more, the correlator then carries continuous parameters, 
namely the positions of the line operators in the extra dimensions. In the context of 
integrable models, continuous parameters on which the Boltzmann weight depends are 
called spectral parameters, and their presence is essential for integrability. Thus, a 2d 
TQFT with extra dimensions produces from line operators a solution of the Yang-Baxter 
equation with spectral parameter. 
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We wish to understand the connection between brane tilings and integrable lattice mod¬ 
els in light of these independent, though apparently related, developments in the relevant 
areas. To this end, it proves helpful to first study the case of the brane box construction [18], 
which is a precursor of the brane tiling construction and conceptually simpler. 

The brane box construction is similar to the M5-brane construction described above. 
In this construction, we start with a stack of N D5-branes on x x T‘^, and put NS5- 
branes on submanifolds of the form x x Ca- Here Ca are 1-cycles of T^, making up a 
lattice. At low energies, this brane configuration realizes a 4d AA = 1 quiver gauge theory, 
which we refer to as a “brane box model,” placed on x 5^. (Actually, can be replaced 
by any Riemann surface without breaking the J\f = 1 supersymmetry.) Just as we did for 
the M5-brane construction, we can relate this theory to a 2d TQFT by considering its 
supersymmetric index. Adapted to the present situation, the argument used there shows 
that the supersymmetric index is given by a correlation function (][][^ Ca{Ca)) in a TQFT 
on r^, where Ca{Ca) is a line operator created by the NS5-brane wrapped around Ca- 

According to the construction of lattice models from TQFTs, this correlator coincides 
with the partition function of a lattice model. To establish the integrability of the model, 
we need extra dimensions along which line operators can move freely. An extra dimension 
indeed emerges as 11th dimension if we embed the brane system into M-theory by string 
dualities. Hence, we conclude that the supersymmetric index of a brane box model is equal 
to the partition function of an integrable lattice model. It turns out that the Yang-Baxter 
equation reduces to Seiberg duality for SU(A^) SQCD with 2N flavors. 

Once the connection between the brane box construction and integrable lattice models 
is understood, the case of brane tilings is not so difficult. In this case we consider defor¬ 
mations of brane box configurations, in which we let the NS5-branes combine with the 
D5-branes over ribbon-shaped neighborhoods of the 1-cycles Ca in T^. Such a deformed 
brane configuration still yields a 4d AA = 1 quiver gauge theory, provided that the defor¬ 
mation meets a certain criterion. By the same reasoning as above, we deduce that the 
supersymmetric index of this theory is given by the correlation function of “thickened” 
line operators representing the ribbon neighborhoods, and equal to the partition function 
of an integrable lattice model. In this way the Bazhanov-Sergeev model arises from brane 
tilings. Again, the Yang-Baxter equation boils down to Seiberg duality, though this time 
the equation involves a sequence of four basic duality transformations. 

So we have answers to the first two of our questions. The supersymmetric indices of 
theories constructed from brane tilings are captured by a lattice model since they are given 
by correlation functions of line operators in a TQFT, which in turn is a consequence of 
the nature of the brane construction and the fact that the index is a protected quantity. 
Furthermore, the integrability of this lattice model is guaranteed by the hidden extra 
dimension which emerges after the brane system is embedded in M-theory. 

In sections 2-4, we discuss in greater detail the connections summarized here among 
the brane box and brane tiling constructions, TQFTs with extra dimensions, and integrable 
lattice models. After reviewing in section 2 the construction of integrable lattice models 
from TQFTs with extra dimensions, we explain in section 3 how it can be applied to the 
supersymmetric indices of brane box models. We identify the associated integrable lattice 


- 3 - 


model, and show that the Yang-Baxter equation for this model takes the form of Seiberg 
duality. In section 4, we treat the case of brane tiling models. 

Sections 5 and 6 are devoted to answering our third question, that is, finding more 
examples of quiver gauge theories whose supersymmetric indices are captured by integrable 
lattice models. 

One way to produce more quiver gauge theories is to apply T-duality to the brane 
configurations discussed above. In this manner we get 3d AA = 2 and 2d M = (2, 2) 
theories. Being related to the 4d parents by T-duality, the supersymmetric indices of these 
theories are also given by the partition functions of some integrable lattice models. In 
section 5, we will see that for these theories, the Yang-Baxter equation follows from lower¬ 
dimensional analogs of Seiberg duality, namely a variant of Aharony duality [19] in the 3d 
case and Hori-Tong duality [20] in the 2d case. 

Perhaps more unexpected is that there are 2d M = (0, 2) quiver gauge theories whose 
supersymmetric indices, or elliptic genera, exhibit integrability. They have half as many su¬ 
percharges as the theories mentioned so far, and cannot be obtained by simple dimensional 
reduction from three or four dimensions. 

In section 6, we discuss three classes of such N = (0, 2) theories. Two of them are 
M = (0, 2) counterparts of the classes oi M = (2, 2) theories considered in section 5. 
The Yang-Baxter equation for these classes identifies two theories related by Seiberg-like 
triality, discovered by Gadde, Gukov and Putrov [21]. The third class consists of A/" = (0, 2) 
theories constructed from brane cube configurations [22], and actually gives rise to a 3d 
lattice model. The integrability condition for this model is that its Boltzmann weight 
satisfies Zamolodchikov’s tetrahedron equation [23, 24], which is the 3d analog of the Yang- 
Baxter equation. Our analysis in this section will be somewhat incomplete, unfortunately. 
For the first two classes, we will demonstrate the integrability of the associated lattice 
models, but not identify the underlying brane constructions or TQFT structures whose 
existence is strongly suggested by the integrability. For the third class, on the other hand, 
we will describe the brane construction and associated lattice model, but not determine 
the corresponding solution of the tetrahedron equation.^ 

Having established the structures of integrable lattice models in several classes of quiver 
gauge theories, we should now ask what we can do with these structures. It is likely that 
knowledge accumulated in the area of integrable models provides new insights into these 
theories or quiver gauge theories in general. Conversely, tools from the gauge theory 
side, such as localization and 1/A expansion, may help further elucidate the physics of 
integrable models. 

Given the generality of the TQFT construction of integrable lattice models, we also 
expect that there are many more applications than those discussed in this paper. Below 
we describe just a few possibilities. 

The NS5-branes in a brane box configuration can be mapped by dualities to either M5- 
or M2-branes intersecting with a stack of M5-branes. These branes represent codimension- 

different solution of the tetrahedron equation has been found recently by Gadde and Yamazaki [25]. 
Their solution is based on A/” = (0, 2) SQCD [21] with all flavor nodes having equal ranks. It remains to be 
seen whether this one may be understood from the perspective adopted in this paper. 
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2 and -4 defects in 6d AA = (2,0) superconformal field theory. From the 4d point of view, 
they are domain walls and line operators in an AA = 2 theory.^ This observation suggests 
that integrability plays a key role in the physics of defects in the 6d theory and 4d A/" = 2 
theories, and the ideas contained in this paper may be useful for understanding it. In this 
regard, we point out that the Yang-Baxter equation has made an unexpected appearance 
in studies of the moduli spaces of AA = 2 theories compactified on [26, 27]. 

Another important appearance of integrability is found in the Bethe/gauge correspon¬ 
dence [28-30] between 2d J\f = (2, 2) theories and quantum integrable systems. It may be 
possible to connect that story to ours, by studying the supersymmetric vacua of AA = (2, 2) 
theories associated with brane box and brane tiling configurations. 

Lastly, integrability features prominently in the AdS/CFT correspondence, and the 
constructions discussed in the present work may shed light on this aspect; after all, the 
developments of brane box and brane tiling techniques were motivated by the AdS/CFT 
correspondence. In fact, brane box models were studied at one-loop level in the planar limit 
in [31], and it was shown that at this level, their dilatation operator can be identified with 
the Hamiltonian of an integrable spin chain. This fact raises the hope that the integrability 
in the AdS/CFT correspondence may be understood in a framework that extends the one 
presented in this paper. 

2 Integrable lattice models from TQFTs with extra dimensions 

As explained in the introduction, the construction of integrable lattice models using line 
operators in TQFTs constitutes an essential ingredient for our argument. So let us begin by 
reviewing this construction. Our discussion mainly follows Costello’s paper [16], to which 
we refer the reader for more details. 

2.1 Vertex models from TQFTs 

Consider a 2d TQFT equipped with a family of line operators parametrized by some 
continuous set, say M or C. We denote by Cr{C) a line operator with parameter r, supported 
on an oriented closed curve C. Let us place this TQFT on a torus T^, and choose 1-cycles 
Ai, ..., Am, Bi, ..., Bn that form an m x n lattice on T^. The case with (m, n) = (2,3) 
is illustrated in figure 1. We are interested in the correlation function 

/ m n \ 

Z{{ra}, = ( n n (Bp)) ( 2 - 1 ) 

'o=l 13=1 ' 

of line operators wrapped around these cycles. 

We compute this correlation function by breaking up the torus into smaller rectangular 
pieces, as shown with dotted lines in figure 1. The idea is that we first perform the path 
integral on each of these pieces, and then glue the results together to reconstruct the path 
integral on T^. 

^We can obtain any M — 2 theory of class S [12, 13] from the brane box configuration (3.1) by replacing 
with the relevant Riemann surface and introducing D5-branes in the 012357 directions. 
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Figure 1: A 2 x 3 lattice of line operators on a torus. The dotted lines divide the lattice 
into rectangular pieces, each of which contains two intersecting segments of line operators. 


Each piece in this decomposition is, topologically, a square with two line operators 
crossing in the middle; 


( 2 . 2 ) 


Using the topological invariance of the theory, we can deform it as 



(2.3) 


Let Vr be the space of states on a finite interval intersected by Cr- Intuitively, Vr is the 
Hilbert space of an open string with a particle attached whose worldline is the line operator 
Cr- In this language, the deformed picture shows an initial state of two open strings in 
Vr ®Vs evolving into a final state in 14 (8> U- that results from interaction in the middle. 
The original picture may be regarded as a particular case where the curved boundaries are 
shrunk to points. Thus, the path integral on the square produces a linear map 


R{r, s): 14 (g) 14 


Vs®Vr. 


(2.4) 


We call R the R-matrix of the TQFT with line operators. 

Choosing a basis {er,i} for 14 for all r, we can represent R by its matrix elements. 
The matrix element (r, s) is the scattering amplitude for a process in which the initial 
state ends up in the final state ^r,i 2 - We represent it pictorially as 



h 

h 

ii 



s 


(2.5) 


- 6 - 









Then, the correlation function (2.1) is given by the formula 


Z{{ra]Asp]) 


m n 


E E nn 

{*a./3} {jp,a} " — 1 /3 —1 


i/3,a+i 

^q;,/3 + 1 

^ck,/3 

i/3,a 


sp 


( 2 . 6 ) 


with the periodic boundary conditions ia,n+i = ia,i and j/ 3 ,m+i = The right-hand 

side is what we get by gluing a collection of pieces like the picture (2.2) to reconstruct the 
line operators wrapped on the torus. 

By representing the R-matrix as above, we are treating the open strings as if their 
physical degrees of freedom are carried solely by the particles attached to them, or equiv¬ 
alently, by the edges of the lattice. In words, the formula (2.6) instructs us to do the 
following. First, choose a state on every edge of the lattice. Next, compute the probabil¬ 
ity amplitude for this configuration of states by multiplying the corresponding R-matrix 
elements. Finally, sum over all possible such configurations to find the answer. 

This procedure is, in fact, precisely how the partition function of a vertex model in 
statistical mechanics is dehned. In a vertex model, state variables are assigned to the 
edges of a lattice, and the interaction takes place at the vertices. The total energy E of the 
system is the sum of interaction energies at the vertices, so the Boltzmann weight 
factorizes into weight factors associated to the vertices. These factors are matrix elements 
of R. In this context, r is called the spectral parameter of the model. 

In conclusion, we have found that the correlation function (2.1) for a lattice of line 
operators in a 2d TQFT coincides with the partition function of a vertex model defined on 
the same lattice, with the line operator parameter playing the role of a spectral parameter. 

2.2 Integrability and extra dimensions 

Let us hx the parameters {s^} in the lattice of line operators under consideration. Viewing 
the vertical direction of the torus as the time direction, we define the Hilbert space of the 
vertex model to be the tensor product^ 

n 

'H = (g)14,. (2.7) 

/3=1 

Furthermore, we introduce the row-to-row transfer matrix 


T{r) = Try^(R(r,si) • • • oy, R{r,Sn)). (2.8) 

®This is the space of states on the disjoint nnion of n intervals, with /3th interval intersected by 
The actual Hilbert space of the TQFT is generally smaller than E, since the degrees of freedom associated 
with the endpoints (“Chan-Paton factors”) must match between adjacent intervals. 
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It is an endomorphism of with matrix elements 


n 


V Jn,l 


(^)=E n 


{ip} /3=1 


jP,2 

ip+i 

ip 

jpp 


sp 


(2.9) 


In terms of the transfer matrix, the partition function (2.6) is written as 

m m 

Z{{ro?i, {s;?}) = Tr^ n = E n (2-10) 

"=1 {jp,a}0‘=^ 

In a TQFT, time evolution is trivial unless the state hits some operator at some point in 
time. We may think of T{r) as the time evolution operator induced by Cr- 

We say that the vertex model on the lattice of line operators is integrahle if T(r) is 
analytic in r and 

[r(r),r(r')] = 0 (2.11) 

for all r, r'. The rationale for this terminology is that by expanding T[r) in powers 
of r, we get a tower of operators that commute with one another and with the time 
evolution operator. 

A sufficient condition for integrability is that the R-matrix has the following two prop¬ 
erties. The first is that R{r, s) is an isomorphism for r ^ s, with the inverse being R{s, r): 

R{r,s)R{s,r) = 1. (2-12) 

The second is that R satishes the Yang-Baxter equation: 

R{s, t)R{r, t)R{r, s) = R{r, s)R{r, t)R{s, t). (2-13) 


Each factor on the two sides of the equation is regarded as an operator on a tensor product 
of Vr, Vs and V). For example, R{r,t) on the left-hand side is a map from Vs®Vr ®Vt 


to Vs®Vt® Vr- 

In the TQFT language, the identity (2.12) means that two tangled line operators can 
be straightened out: 



(2.14) 


r 


s 


r 


s 
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The Yang-Baxter equation, on the other hand, expresses invariance under moving a line 
operator past the intersection of two other line operators: 



r s 


r s t 



(2.15) 


To see that these properties imply the commutativity of transfer matrices, consider the 
following move of line operators: 



(2.16) 


Taking the trace over Vr ® V).' of both sides, and using the identity (2.14) and the cyclic 
property of the trace, we find that the left- and right-hand sides become T{r)T{r') and 
T(r')T(r), respectively. 

In general, a 2d TQFT does not have the above properties. Physics may change 
discontinuously under the relevant moves due to phase transitions, as the topology of 
the line operators does not remain the same. Imagine, however, that there are “extra 
dimensions” hidden in the above pictures. For instance, the direction perpendicular to 
the page may provide one. The line operators can then sit at different points in the 
extra dimensions. If so, these moves are topologically trivial and hence the R-matrix 
satisfies the desired properties. An argument of this sort is familiar from studies of knot 
invariants, which may be approached either from the perspective of 3d TQFTs or lattice 
models [32, 33]. 

The above observation motivates us to introduce the notion of extra dimensions to 
TQFTs. Consider a L>-dimensional quantum field theory formulated on S x M, where 
S is any d-manifold and M is a fixed [D — (i)-manifold. We say that the theory is a 
d-dimensional TQFT with extra dimensions if it is topological on We refer to M as 
the internal space of the theory. This definition is very much in the same spirit as the 
notion of a 2d conformal field theory “valued in 4d quantum field theories” [47] developed 
in connection to the 6d construction of 4d A/" = 2 theories of class S [12, 13]. In that case, 

^Although such theories might sound exotic, actually many have already been studied before. Typically, 
they are constructed from supersymmetric theories by topological twisting along S. 2d examples are 
Costello’s theory mentioned below, 4d A/” = 2 superconformal gauge theories with M being any Riemann 
surface [34], 5d A/” = 2 super Yang-Mills theory with M = [35], and 6d A/” = (2, 0) superconformal field 

theory with M = x [14]. 3d examples include 5d A/” = 2 super Yang-Mills theory with M = [36, 37] 

and [38-40], 6d A/” = (2, 0) theory with M = S^xS^ [41], [42], and more general lens spaces L(j), q) [43]. 

4d examples are 6d A/” = (2,0) theory with M any Riemann surface [44-46]. 
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however, the choice of the 4d spacetime is not part of the definition. By contrast, here we 
demand the theory to be topologically invariant on S for a specific choice of M. 

Now we come to the main point of our discussion. Suppose that a 2d TQFT with extra 
dimensions has line operators. Place the theory on and wrap line operators around 1- 
cycles Aa and Bp. In the internal space M, the 1-cycles are located at some points Xa and 
yp, respectively. The situation is similar to what we considered at the beginning of this 
section, but unlike that case, this time we do not assume that the line operators carry a 
continuous parameter. To highlight the difference, let us consider the extreme case that 
the theory has only one type of line operator, C. Repeating the same argument, we arrive 
at the conclusion that the correlation function 

/ m n \ 

z {{ xa } dy ^]) = (n C{Aa X { Xa }) n£(B^x{M))) (2.17) 

'a=l p=l ' 

is given by the partition function of a vertex model defined on the lattice of line operators. 
As we have just seen, the presence of extra dimensions guarantees that the R-matrix is 
invertible and satisfies the Yang-Baxter equation. 

A beautiful insight of Costello [16, 17] is that even though the line operator lacks a 
continuous parameter, the R-matrix still depends on such parameters: the positions Xq, yp 
of line operators in M. Thus, the vertex model is integrable if the transfer matrix varies 
analytically on M. Based on this idea, it was shown in [17] that a special case of the 
6-vertex model corresponding to the XXX spin chain, and its generalizations based on Lie 
algebras other than 5 I 2 , arise from a 2d TQFT with extra dimensions whose internal space 
M = CP^. This TQFT is obtained from a deformed and topologically twisted version of 4d 
M = 1 super Yang-Mills theory, and the line operators used in the construction are Wilson 
lines. The extra dimensions therefore elegantly explain not only why the Yang-Baxter 
equation holds for these models, but also where the spectral parameter comes from. 

Speaking of extra dimensions, it should be noted that although only line operators 
were considered in our discussion, we could as well use higher-dimensional operators that 
have codimension at least two in S x M. This is because after we wrap them on 1-cycles 
on S, there is still room in M for them to avoid one another. In our main examples, we 
will make use of codimension-2 defects in a 6d theory. 

2.3 IRF models 

In the construction described above, we obtained a vertex model from line operators by 
letting the midpoints of open strings represent all physical degrees of freedom. In many 
cases, however, degrees of freedom really reside on (and only on) the endpoints. In such 
cases it is more efficient to reformulate the vertex model as an interaction-round-a-face 
(IRF) model. 

In an IRF model, state variables live on the vertices of a lattice, and interaction takes 
place among vertices connected by edges that surround a face. We use letters a, b, ... to 
denote state variables. For a square lattice, the Boltzmann weight for the interaction is 


10 


denoted as 


W 


/ d c 
\a b 


r,s = 


©- 




■ 4 ) 


(2.18) 


The dashed oriented lines are rapidity lines, and make up the dual lattice. The parameters 
r, s are spectral parameters, also called rapidities. An IRF model is a vertex model on 
the dual lattice, in which the state space assigned to an edge is the direct product of two 
spaces and many elements of the R-matrix vanish. 

From the TQFT point of view, rapidity lines are line operators, and the rapidities are 
their parameters. The Yang-Baxter equation (2.15) for an IRF model is 




(2.19) 


In terms of the Boltzmann weight, the equation reads 


9 


e d 
9 c 


s,t]W 


9 c 
a b 


r,t]W 


e 9 
f a 




r, s 


W 


e d 
f 9 


r,t]W 


f 9 
a b 



( 2 . 20 ) 


2.4 3d TQFTs and the tetrahedron equation 


As a generalization, we can construct higher-dimensional lattice models using codimension- 
1 defects in a d-dimensional TQFT with d > 2. Placing the theory on a d-torus and 
putting these defects on (d — l)-cycles that form a d-dimensional lattice, we get a vertex 
model on this lattice. If the theory has extra dimensions, the R-matrix of the model satisfies 
a d-dimensional analog of the Yang-Baxter equation. 

To be concrete, let us take d = 3. In this case, codimension-1 defects are surface 
operators, and three of them intersect at a vertex of the lattice. A cubic neighborhood of 
a vertex looks like 



( 2 . 21 ) 
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The orientations of the lattice edges specify those of the surface operators (drawn with 
dotted lines in the picture). Viewing this picture as the worldvolumes of three faces of 
the cube traveling in a diagonal direction, we see that the R-matrix R{r, s, t) produced by 
the path integral on this cube is an endomorphism of V-,s ® 14 ,t <8) 14 ,^, where 14 ,* is the 
Hilbert spaces on a square intersected by two surface operators with parameters r and s. 
Pictorially, the R-matrix is represented as 


R{r,s,t) 



( 2 . 22 ) 


To avoid clutter, in what follows we will suppress spectral parameters in this kind 
of pictures. 

The 3d analog of the Yang-Baxter equation is Zamolodchikov’s tetrahedron equa¬ 
tion [23, 24], obtained by moving one of four surface operators that form a tetrahedron: 



(2.23) 


The tetrahedron equation implies the commutativity of layer-to-layer transfer matrices at 
different values of the spectral parameter. 

When the state variables may be considered as living on the cubes surrounded by 
surface operators, a vertex model can be reformulated as an interaction-round-a-cube (IRC) 
model [48]. The Boltzmann weight of an IRC model is denoted as 



The tetrahedron equation for an IRC model is 


^ lT(a 4 lci, C3, C2I63, b2,bi\d)W{ci\a3, 61,62]^, cq, Ci\bi) 
d 

X IT(6i1c4, C3, d\b3, 64, a2lc5)lT((il64,63, b 2 \c 2 , cq, csjoi) 

= ^ 1T(6i1c 4, C3, cila4, 03, a2ld)lV(cila3, 04, 62IC2, cq, d\ai) 

d (2.25) 

X W{ai\d, C3, C2I63, ai, a2lc5)lT(^[03, 02, oijcs, cq, C4I64). 


Graphically, it takes the form of equivalence between two ways of dividing a dodecahedron 
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into four hexahedra: 




(2.26) 


3 4d A/” = 1 quiver gauge theories 

Now we wish to apply the above construction to 4d AA = 1 quiver gauge theories realized 
by brane box configurations. In this section, we explain how the elements that entered the 
construction — TQFT, lattice of line operators, and extra dimensions — arise nicely in the 
supersymmetric indices of brane box models. Furthermore, we determine the associated 
integrable lattice model, and demonstrate that the Yang-Baxter equation for this model 
is equivalent to the equality between the supersymmetric index of S\J{N) SQCD with 2N 
flavors and that of its Seiberg dual. 

3.1 Brane box models 

Consider Type IIB superstring theory, and suppose that we have N coincident D5-branes 
extending in the directions 012346, and a number of NS5-branes in 012345 and 012367: 
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These branes are all located at the same point in the 89-plane. So each NS5-brane intersects 
with the D5-branes along a straight line parallel to the 4- or 6-axis. 

Macroscopically, this brane system realizes a supersymmetric quiver gauge theory, 
much the same way as the brane constructions of 3d AA = 4 [49] and 4d AA = 2 [11] gauge 
theories. In those constructions, D3- or D4-branes are suspended between NS5-branes 
placed at finite intervals in the 6-direction. The difference here is that the D5-branes are 
divided by NS5-branes into finite segments in the 4-direction as well. The resulting theory 
therefore has less spacetime dimensions and supersymmetry, but a more complicated quiver: 
it is a 4d AA = 1 gauge theory described by a planar quiver diagram. Brane box models, 
first studied in [18], are theories constructed from brane configurations of this kind. 

To avoid complications associated with semi-infinite branes, let us compactify the 46- 
plane to a torus T^. Then, the intersections of the NS5-branes with the D5-branes form 
a lattice on the 46-torus. We draw the intersections by oriented dashed lines and call 
them rapidity lines, anticipating connection with integrable lattice models. We choose the 
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Figure 2: The 2x3 brane box model, (a) The configuration of NS5-brane intersections 
on the 46-torus, (b) The corresponding quiver diagram. 


orientations for horizontal lines in such a way that they are parallel (and not antiparallel) 
to one another, and similarly for vertical lines. See figure 2(a) for an example of a brane 
box model, in which the rapidity lines make up a 2 x 3 lattice. 

The rule for determining the quiver is as follows. We label the faces of the lattice by 
letters a, b, .... To the ath face is associated a vector multiplet for gauge group SU(iV)a, 
which originates from open strings whose ends are attached on this part of the D5-branes. 
Since the rank of the gauge group is fixed throughout our discussion, in our quiver diagram 
we represent the vector multiplet simply by a circle node labeled a: 

@ (3.2) 

From two faces a and b adjacent to each other or sharing a vertex, we get a chiral multiplet 
in a bifundamental representation of SU(iV)a x SU(A^);,. This is represented by an arrow 
between the two nodes: 

(a)-(3.3) 

is a chiral multiplet in the representation (N, N) of SU(A^)a x S\J{N)h, coupled to the 
vector multiplets. If four faces a, b, c, d are placed around a vertex in the counterclockwise 
order, the corresponding part of the quiver is given by 



What makes brane box models physically interesting is that they are chiral: between nodes, 
arrows only point in one direction. The quiver diagram for the 2x3 brane box model is 
shown in figure 2(b). 

The same diagram also encodes the superpotential of the theory. The quiver consists of 
a collection of triangles that cover the 46-torus. For each triangle formed by three arrows, 
there is a cubic superpotential term given by the trace of the product of the corresponding 
bifundamental multiplets, with sign determined by the orientation of the arrows. The 
precise form of the superpotential will not be important for our discussion; the quantities 
we will consider are largely independent of this information. 
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Figure 3: Sliding D5-brane segments along NS5-branes. Illustrated here is an ideal situ¬ 
ation with ffs = 0. 


NS5 NS5 



(a) (b) 


Figure 4: For gg > 0, taking apart two D5-brane segments creates either (a) a tilted {N, 1) 
5-brane or (b) a tilted {N, —1) 5-brane. 

The gauge group of a brane box model is the product of SU(A^)s and not U(A^)s. To 
understand why, let us suppose for simplicity that there are only NS5-branes spanning the 
012345 directions. In this case we get a 5d M = 1 theory. We are taking the string coupling 
Qs to be small, which means that the NS5-branes have tension much larger than that of 
the D5-branes and can be treated as rigid objects. These heavy NS5-branes “chop” the 
D5-branes into finite segments in the 6-direction. Imagine sliding these segments slightly in 
the 5-direction. If gg were strictly zero, the NS5-branes would be infinitely heavier than the 
D5-branes and not affected by such a deformation at all. Then we would have a situation 
sketched in figure 3. 

In reality, gg is small but nonzero, and the deformation does affect the NS5-branes. 
What actually happens is that when a stack of N D5-branes hit on an NS5-brane, these 
branes combine to form a bound state — an (iV, 1) or {N, —1) 5-brane, depending on the 
relative positions of the branes — which then goes off diagonally, as shown in figure 4. As a 
result, moving the D5-brane segments in the 5-direction entails shifts of semi-infinite NS5- 
branes in the 6-direction, and this costs an inhnite amount of energy. The x^-coordinates 
of the D5-brane segments are therefore frozen. As the diagonal U(l) gauge field on each 
segment is in the same 5d multiplet as the scalar x^, it is frozen as well. Thus, the gauge 
group of the 5d theory is the product of SU(A^)s, times a U(l) factor corresponding to the 
overall center of mass coordinate of the D5-branes. The U(l) factor is decoupled from the 
rest of the dynamics, so we ignore it. 
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3.2 R-symmetry and flavor symmetries 

The brane system (3.1) is invariant under rotations in the 89-plane. The rotational sym¬ 
metry appears in the effective 4d theory as an R-symmetry U(l)/j. However, the R-charges 
of the bifundmanetal multiplets are not determined uniquely by the brane conhguration, 
because of the possibility of shifting them by U(l) flavor charges. This ambiguity plays a 
crucial role in the connection to integrable lattice models, so let us look at it closely. 

We denote by Rab the R-charge of the bifundamental multiplet represented by an arrow 
going from node o to b. There are two conditions for a set of R-charges {Rab} to be realized 
in a brane box model. The first is that the superpotential must have R-charge 2 so that it 
preserves U(1 )r. Thus, for every triangle 



we must have 

Rab + Rbc + Rea = 2. (3.6) 

The second is that U(l)ij must be nonanomalous. This requirement leads to the condition 
that at every node, the R-charges of the bifundamental multiplets starting from or ending 
on that node add up to the number of such multiplets minus two. In our case, for every 



we must have 

Rab + Rea + Rad + Rea + Raf + Rga = 4. (3.8) 

The second condition is equivalent to the vanishing of the one-loop /^-functions for the 
gauge couplings, and one expects that the theory flows to a nontrivial infrared fixed point 
when this condition is satisfied. 

What is the dimension of the space of solutions to the constraint equations? The 
quiver of a brane box model provides a triangulation of T^. The number of R-charges to 
be assigned is equal to the number of edges E, whereas the number of constraints is equal 
to the number of faces F plus that of vertices V. Thus, naive counting would suggest that 
the dimension is given by minus the Euler characteristic y = R — E -|- E, which vanishes 
for T^. This is not the case, though. For example, we can assign R-charge i? = 1 — r to 
the horizontal arrows, E = 1 -|- r to the vertical ones and E = 0 to the diagonal ones, with 
r being a free parameter. 

The reason for the mismatch is that the constraints are not all independent. In¬ 
deed, using the triangle condition (3.6), the anomaly cancellation condition (3.8) can be 
replaced by 

Rad ~ R-be + Rga ~ Rfe = 0 - ( 3 - 9 ) 
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Figure 5: The R-charge parametrization for brane box models. 


The sum of the left-hand side, with a running over all nodes in the same row of the quiver, 
is zero due to the periodicity of . So there is one relation for each row. Similarly, there is 
one relation for each column. For an m x n lattice, the dimension of the space of solutions 
is therefore at least m + n, the total number of rapidity lines. 

These m + n degrees of freedom in the R-charge assignment may be thought of as 
degrees of freedom associated to the rapidity lines. Concretely, we can parametrize the 
R-charges as follows. To each rapidity line we assign a real parameter, and call it the 
rapidity or spectral parameter of the rapidity line. Then, an arrow has R = 1 — r if it is 
horizontal and R = 1 -|- r if vertical, where r is the spectral parameter of the rapidity line 
intersecting with that arrow. The R-charges of the diagonal arrows are determined by the 
condition (3.6). The rule is summarized in figure 5. One readily sees that with this rule, 
the condition (3.8) is also satished. 

In general, a U(l) R-symmetry can be shifted by U(l) flavor symmetries, that is, global 
U(l) symmetries that commute with the supersymmetry algebra. In view of this fact, the 
degrees of freedom in the R-charge assignment just found suggests that there are m -|- n 
U(l) flavor symmetries, each associated to a single rapidity line. 

The origin of these global symmetries can be most clearly understood in the original 
brane picture. Let us pick an NS5-brane, say the one corresponding to the ath rapidity 
line, and imagine, as before, slightly sliding the D5-brane segments that end on it from left 
and right. This operation splits the NS5-brane into the upper and lower semi-inhnite NS5- 
branes, joined by a finite segment of {N, ±1) 5-brane; see hgure 4. For definiteness, consider 
the case with an {N, 1) 5-brane. Bifundamental multiplets arise from open strings whose 
ends are localized at the junctions where semi-infinite NS5-branes meet D5-branes and the 
{N, 1) 5-brane. These strings may be viewed as fundamental strings ((A^, 0)-strings), as in 
figure 6(a), or D-strings ((0, l)-strings) with the opposite orientation,^ as in figure 6(b). As 
such, they are acted upon by global U(l) gauge transformations on the semi-infinite NS5- 
branes. In the effective 4d theory that emerges when the split branes are rejoined, gauge 
transformations on the upper and lower branes appear as elements in the same symmetry 

®To determine the orientation, we draw an {N, 0)-string attached on the D5-branes, and move the end¬ 
points to the semi-infinite NS5-branes. When an endpoint passes a junction, an {N, l)-string is created due 
to charge conservation. When the other endpoint passes the other junction, the {N, l)-string is annihilated, 
leaving a (0, l)-string. The orientation of the (0, l)-string can be seen from the final picture. 
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Figure 6: (a) An (A^, 0)-string stretched between the D5-NS5-(A^, 1) junctions, (b) The 
same string may be regarded as a (0, l)-string with the opposite orientation. 


U(l)o, just acting in opposite manners. Hence, a U(l) flavor symmetry is associated to 
each NS5-brane or rapidity line. 

Equivalently, we may regard U(1 )q, as arising from the U(l) gauge symmetries on the 
D5-brane segments, which we have seen are nondynamical. From this point of view, U(1 )q, 
acts on all D5-brane segments on the “right” side of the ath N5-brane, or on all of those 
on the “left” side in the opposite manner. Of course, the words “left” and “right” are 
not really meaningful when the 4- and 6-direction are periodic, but we can make sense 
of this interpretation by, roughly speaking, considering the D5-branes to be sections of a 
U(l)Q,-bundle on the 46-torus. 

This point is probably easier to understand if we consider the situation where there are 
only NS5-branes in the 012345 directions, and look at the mass parameter associated 
to U(l)a in the corresponding 5d AA = 1 theory, nia is the value of the real scalar in the 
background vector multiplet for U(1 )q, and proportional to the difference between the x^- 
coordinates of the two D5-brane segments that end on the ath NS5-brane. The periodicity 
in the 6-direction suggests that the sum of over the rapidity lines must vanish. However, 
this constraint can be modified if we compactify the 6-direction using a twisted boundary 
condition that involves a shift in the 5-direction. (The same twisted boundary condition 
is used in the brane construction of 4d AA = 2 elliptic models [11].) This modification 
introduces one more parameter, namely the amount of shift, leading to the same number 
of mass parameters as the rapidity lines. 

We normalize the U(l)Q-charge in such a way that horizontal and vertical bifunda¬ 
mental arrows intersecting with the ath rapidity line have = —1 and -t-1, respectively, 
while Fjj = 0 for {3 ^ a. The charges of the diagonal arrows are found by demanding that 
the superpotential have F^, = 0 for all a. Using F^, our R-charge parametrization can be 
concisely written as 

R = Ro + '^ VaFa, (3.10) 

a 

where Rq is the R-charge defined by assigning Rq = 1 to the horizontal and vertical arrows, 
and Tq, is the parameter of the ath rapidity line. As expected, the degrees of freedom in 
the definition of R come from shifts by flavor symmetries. 

We have seen above that there are at least m + n U(l) flavor symmetries, but have not 
determined the precise number. Actually, in some cases there are extra U(l) symmetries in 


- 18 - 









Figure 7 : The general R-charge assignment for the 2x2 brane box model. Here only the 
R-charges of the horizontal and vertical arrows are indicated, and the variables satisfy the 
relation r' — r = s — s' = t — t' = u' — u. 

addition to the ones discussed above. For example, for a 2 x 2 lattice, the most general R- 
charge assignment has five parameters rather than four, as shown in hgure 7. In fact, from 
the relation between the brane box and brane tiling constructions explained in section 4, 
one can deduce that the number of U(l) flavor symmetries is equal to m+n+gcd(m, n) — 1. 
We will not consider these extra flavor symmetries. Alternatively, one may assume that m 
and n are coprime. 

3.3 Structure of a TQFT with extra dimension 

So far we have described the brane box configuration (3.1) from the point of view of the 
effective 4d AA = 1 theory. To make contact with the framework discussed in section 2, 
let us instead think of it as realizing 6d maximally supersymmetric Yang-Mills theory on 
the D5-branes, in the presence of codimension-1 defects created by the NS5-branes. We 
wish to turn this 6d theory into a theory that is topological on the 46-torus, and regard 
the defects as line operators with spectral parameter in the TQFT. This is achieved by 
applying an appropriate topological twist and restricting physical quantities to compute, 
as we now explain. 

Type IIB superstring theory has M = (2, 0) supersymmetry in 9 -|- 1 dimensions, 
generated by Majorana-Weyl spinors ei, eR such that 

rol23456789eL = roi23456789ei? = ^R- (3-11) 

Here Tm, M = 1, ..., 9 are gamma matrices, and roi 23456789 = Tq • • • Tg is the chirality 
operator. The D5-branes preserve supersymmetry generated hy er, er satisfying 



ei, = roi2346ei?- 

(3.12) 

Likewise, 

the NS5-branes impose conditions 




= roi2345eL, 

ER — —roi2345eR 

(3.13) 

and 

Cl = roi2367eL, 

£« = —roi2367eR 

(3.14) 
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on the preserved supersymmetry. The conditions on coming from the NS5-branes ac¬ 
tually follow from those on and the D5-brane condition, hence are redundant. Since 
roi2345 and roi2367 Square to the identity and are traceless, half of their eigenvalues are -|-1 
and the other half are —1. Moreover, they commute with each other and are simultaneously 
diagonalizable. Each NS5-brane condition therefore halves the number of independent so¬ 
lutions for e^. Given e/,, the D5-brane condition determines e/j. In total, the 5-branes 
preserve four supercharges, which generate the M = 1 supersymmetry of the 4d theory we 
have been studying. 

From the NS5-brane conditions, we see or 

(r46-I-r57)eL = 0. (3.15) 

This equations says that if we replace the rotation group 80(2)46 of the 46-plane (which 
we decompactify for a moment) with the diagonal subgroup 80 ( 2)45 of 80 ( 2)46 x 80 ( 2 ) 57 , 
then the preserved supercharges become scalars under the new rotation group. From the 
same equation we also find that the N85-brane conditions on el can be generalized to 

CL = roi23(r4Cos0 -I- r6sin0)(r5Cos0 -I- r7sin0)eL, (3.16) 

where 6 is any angle. This shows that the N85-branes can be rotated in the 46- and 57- 
plane by the same angle 6 without destroying the supercharges, or for that matter, they can 
intersect with the 46-plane along arbitrary closed curves as long as the slopes are correlated 
between the 46- and 57-plane. 

From the point of view of the 6 d theory on the D5-branes, the replacement of 80(2)46 
with 80 ( 2)45 is a topological twist along the 46-plane. After the twisting, eight out of the 
sixteen supercharges are scalars on the 46-plane, and as such, can be preserved even when 
the 46-plane is replaced by a general 2-manifold S. Codimension-1 defects can be inserted 
along arbitrary closed curves on S while preserving four of the eight supercharges. These 
defects carry a spectral parameter, as we have seen from the 4d point of view. 

Now, to produce a desired TQFT from this twisted theory, we simply restrict ourselves 
to quantities that are independent of the metric on S, and depend on the codimension -1 
defects inserted on S only through their topology and spectral parameters. Then, the 
twisted theory becomes a TQFT on S, and the defects become line operators with spectral 
parameter in this TQFT. 

A good example of such a topological quantity is the supersymmetric index of the 
corresponding 4d AA = 1 theory [50-52]. As the theory has the R-symmetry U(l)i 7 , it can 
be placed on the Euclidean spacetime x without breaking supersymmetry [50, 52, 53]. 
Under (the double cover of) the isometry group 8U(2) x 8U(2) of the four supercharges 
transform as a pair of doublets (2,1). Let Ji, J' be the generators of the 8U(2) factors. 
The supersymmetric index is dehned by 

X{p,q,M,{ra}) = TV53('(-1)V^+^3+«/2/3-J^+R/2 

where p, ( 7 , Ua are complex parameters, and the trace is taken over the Hilbert space on 
S^. This quantity is topological for the following reason. Bosonic and fermionic states 
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with energy E ^ 2 J 3 + K are paired by the action of a supercharge. (Here we are taking 
the radius of to be 1.) Due to the presence of the factor (—1)^, pairwise cancellations 
occur and the only contributions to the trace come from states with E = 2 J 3 + E. Under 
continuous changes of the parameters of the theory (other than the rapidities on which 
K implicitly depends), the index remains invariant since states can be brought into or out 
of the right energy level only in boson-fermion pairs. The index is therefore protected 
against such changes, in particular against deformations of geometric parameters such as 
the shapes of line operators and the metric on S. 

There are other topological quantities. For example, we can replace the in the above 
definition of the supersymmeric index with x or a lens space L{p,q) [54]. Unfortu¬ 
nately, though, an explicit expression for the index is not known except for simple cases 
such as L{p, 1) [15]. We will focus on the index, for which the relevant mathematical 
results are available. 

So we have obtained from brane box models a 2d TQFT that has a family of line op¬ 
erators with spectral parameter. We can now apply the construction discussed in section 2 
and get a vertex model, by placing the TQFT on and making a lattice of line operators. 
However, we cannot tell whether the model is integrable just from the structure of a TQFT. 
For the integrability, it is necessary to show, in addition, that no phase transitions occur 
when two line operator are untangled or a line operator passes through the intersection of 
two other line operators. In Costello’s argument, these properties are guaranteed by the 
existence of extra dimensions in which the line operators can miss one another. 

We can make a similar argument in the present case. Let us go back to the brane 
picture. To the configuration (3.1), we apply T-duality in the 3-direction (which we take 
to be the S'^-direction of the spacetime x 5^) to convert the D5-branes to D4-branes, 
and then lift the resulting Type HA brane system to M-theory: 
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These dualities transform the 6 d theory on the D5-branes into an AA = (2,0) theory on 
a stack of M5-branes, and the codimension-1 defects into codimension-2 defects. There is 
now an extra dimension along which line operators can move, namely the M-theory circle. 
Its presence ensures the integrability of the vertex model. 

The similarity to Costello’s argument suggests that the x^*^-coordinates of line op¬ 
erators are related to the spectral parameters r^. This is indeed true, as the following 
analysis shows. 

If we trade off the flavor fugacities Ua for new parameters Sq given by Ua = (m)®“^^, the 
supersymmetric index (3.17) depends on Sq, only through the combination R + SaEa- 
Hence, the effect of turning on nonzero values for Sa is equivalent to shifting Tq, —>■ Tq, -|- Sq, 
and Sa represent the same degrees of freedom as Tq,. We want to understand how Sa are 
described in the language of branes, and keep track of them under the dualities. 
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As we saw before by considering the deformation depicted in figure 6, U(l)a can be 
interpreted as either the global gauge symmetry on the lower half of the ath NS5-brane, or 
that on the upper half acting in the opposite manner. The conserved current Ja for U(1 )q 
is therefore supported along Ca on the 46-torus, and couples to the 6d theory through a 
term of the form 


i I (^Q ^q) 

Js^xS^xT^ 


(3.19) 


Here A~ and A+ are the gauge fields on the lower and upper halves of the NS5-brane, 
restricted to the brane intersection; they are related there to the U(l) parts of the gauge 
fields on the D5-brane segments by some boundary condition. We give a background value 
to the gauge fields such that A~ — A+ is constant along and Ca- As Fa is the integral 
of -kJa over and T^, we find that the above term becomes 


i 



{A--At) Fa- 


(3.20) 


This term enters the path integral as an exponentiated factor, so Sq, is proportional to the 
difference between the holonomies of A~ and At around the 3-circle. 

These holonomies are mapped under the T-duality to the background values of a scalar 
field on the upper and lower halves of the corresponding Type IIA NS5-brane. Upon lifting 
to M-theory, the NS5-brane becomes an M5-brane, and the scalar is mapped to the x^^- 
coordinate. Therefore, the spectral parameter for a line operator is given by the difference 
between the x^^^-coordinates of the upper and lower halves of the M5-brane. 

We remark that if the holonomies of At and A~ are equal, then the whole ath NS5- 
brane is mapped to a single point on the M-theory circle, and the supersymmetric index 
does not depend on its position at all. This observation is consistent with the fact that 
the J\f = (2,0) theory compactihed on has a topological sector that is equivalent to 
Chern-Simons theory with complex gauge group at level k = 1 [40]. In the present setting, 
the NS5-brane creates a Wilson line in the Chern-Simons theory on x 5^, and its precise 
position on the 3-manifold is irrelevant. When the holonomies are different, the Wilson 
line splits into two line operators, and apparently they no longer belong to the topological 
sector. It is clear from our viewpoint that the same consideration should apply to other 
indices obtained by replacing with appropriate 3-manifolds. Indeed, the M = (2,0) 
theory on x has a topological sector equivalent to complex Chern-Simons theory at 
level k = 0 [38, 39], while the lens space L{k, 1) gives complex Chern-Simons theory at level 
k [43] . Similarly, a twisted product of and leads to analytically continued Chern- 
Simons theory, which is the holomorphic sector of complex Chern-Simons theory [37, 55-57]. 


3.4 Integrable lattice model for brane box models 

Finally, let us identify the integrable lattice model associated with brane box models. To 
do so, we look at the integral formula for the supersymmetric index. 

Recall that the flavor fugacities Ua represent the same degrees of freedom as the spectral 
parameters r^. Since we are going to keep the dependence on r^, we set = 1. With this 
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understood, the supersymmetric index of a brane box model is computed as follows [50, 58]. 
First, we assign the factor 




{{p-,p)oo{q]q)oo)^ 




T{Za,i/Za,j-,P, q) 


to the vector multiplet represented by node a, and 


X^{Za,Zh-,Rab) = X\'^{{pq)^‘''’^‘^Zh^3/Za,i]P,q) 

ij 


(3.21) 


(3.22) 


to the bifundamental chiral multiplet with R = R^b represented by an arrow from node a 
to b. Here 


T{z-,p,q) 


n 


i,j=0 


I — Z 

1 — zp^q^ 


(3.23) 


is the elliptic gamma function, (z; g)oo = n“o(i — zq^), and Za collectively denotes the 
fugacities Za,i, i = 1, ..., N ioi the gauge group SU(A^)a, obeying ^a.,i = 1- Then, we 
multiply the factors from all nodes and arrows, and integrate over the fugacities: 


1 = 


n^n 


Af-l 




J =1 


2 TTiz, 


a, 3 


Xv^^a) X]^(^Za, Zb, Rab) • 

@—*'® 


(3.24) 


The integration for each fugacity is performed over the unit circle = 1}. 

The above procedure clearly shows the statistical mechanical nature of the super- 
symmetric index. In the statistical mechanical interpretation, Za,i are continuous state 
variables, and the factors Xy and Xb are the Boltzmann weights for self- and nearest- 
neighbor interactions, respectively. The partition function is defined by multiplying the 
weights and summing over all possible state configurations, and this is identified with the 
supersymmetric index. 

Since the state variables Za,i are assigned to the nodes of the quiver which correspond 
to the faces of the lattice on the 46-torus, it is natural to formulate this model as an IRF 
model. We can take its Boltzmann weight to be® 


W 


I Zd Zc 
\Za Zb 


r,S j = Xy(Za)XB(Za, Zb] 1 - s)XB{zb, Zc] 1 + r)XB(^c, Za] S - r). (3.26) 


°This is not the only choice. For example, one may replace the right-hand side with the more symmetric 
expression 

\/lw{Za)XB{Za, Zf, 1 “ s)lB{Zb, Zc, 1 "t" r)lB{Zc, Za\ S - v) 

X yiv(2c)XB(2c,2a; s - r)lB{za,Zd-, 1 -|- r)lB{Zd, zy, 1 - s) (3.25) 


without spoiling the integrability. 
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Graphically, it can be represented as 





S 


(3.27) 


(Notice a slight discrepancy in the notation compared to the original definition (2.18). Here 
a, b, etc. on the left-hand side really stands for Za, Zf,, etc.) A square node indicates 
that the associated factor Iy{za) is absent, and a node (d) drawn with a dashed line means 
that the factor is present but integration is not performed over Za- The intersections of 
NS5-branes with the D5-branes, which create line operators in the TQFT on the 46-torus, 
play the roles of rapidity lines, as anticipated. 

Let us show that this Boltzmann weight satisfies the Yang-Baxter equation (2.19). For 
this purpose we will need the identity 


TB{Za, Zb] Rab)d^B{zb, Za] Rba) = 1 (3.28) 

which holds for Rab + Rba = 2. In effect, it means that we can erase two arrows going in 
opposite directions between the same pair of nodes: 



(3.29) 


This property reflects the fact that if one has chiral multipets ^hi, $2 in conjugate rep¬ 
resentations whose R-charges add up to 2, one can give them masses by turning on the 
superpotential term The supersymmetric index is independent of m, and in the 

limit m —>■ 00 , these multiplets become inhnitely massive and decouple, leaving no contri¬ 
bution to the index. 

Now, plugging the expression (3.27) into the Yang-Baxter equation (2.19), we find 




Some nodes and arrows cancel out between the two sides, while the arrows between nodes 
h and g on the right-hand side cancel thanks to the identity (3.29). Moving the arrows 
that are not connected to node g from the left- to right-hand side using the same identity, 


- 24 - 



















we end up with the equality 



(3.31) 


Interestingly, this equality expresses precisely the invariance of the supersymmetric index 
under Seiberg duality. 

Seiberg duality says that for Nf > Nc + 1, an Af = 1 SU(iVc) gauge theory with Nf 
flavors (Q, Q) and no superpotential has a dual description in the infrared as an AA = 1 
SU(A^/ — Nc) gauge theory with Nf flavors {q,q) and Nj singlets M, coupled through a 
superpotential W = qMq. For both theories, the flavor symmetry is SU(A^/)l x S\J{Nf)ji x 
U(l)s and the R-symmetry is U(l)ij. Under these symmetries and the gauge group G = 
SU(A^c) or S\J{Nf — Nc), the chiral multiplets transform as follows: 



G 

SU(iV/)L 

S\J{Nf)R 

U(1)b 

U(1)r 

Q 

□ 

□ 

1 

1 

1 - Nc/Nf 

Q 

□ 

1 

□ 

-1 

1 - Nc/Nf 

q 

□ 

□ 

1 

-Nc/iNf - Nc) 

Nc/Nf 

q 

□ 

1 

□ 

Nc/{Nf - Nc) 

Nc/Nf 

M 

1 

□ 

□ 

0 

2(1 - Nc/Nf) 


(3.32) 


In the case at hand, Nf = 2Nc = 2N and the charged chiral multiplets all have R = 1/2. 
This corresponds to the situation with (r, s,t) = (—1/2, 0,1/2). As we stressed, however, we 
can shift the R-charge by U(l) flavor charges. The quivers that appear in the equality (3.31) 
split the 2N flavors into a pair of N flavors. This breaks each of S\J{2N)l and S\J{2N)^ 
down to an SU(A') x SU(A') x U(l) subgroup. So in total there are three U(l) flavor 
symmetries that can be used to shift the R-charge, and they give rise to the three parameters 
r, s and t. 

Since the supersymmetric index is protected from continuous changes of the parameters 
of the theory, it is invariant under renormalization group flow and therefore under Seiberg 
duality. Indeed, as observed in [58], the invariance under Seiberg duality is equivalent to 
an identity for the elliptic gamma function proved in [59].^ Thus, we have demonstrated 
the integrability of our model.® 


^Historically, the Identity was first discovered for Nc = 2 and Nf = 3 in [60]. Subsequently it was proved 
for Nc = 2 and Nf = 4 in [61] and extended to the general case In [59]. 


“The invertibility of the R-matrix follows from the identity 


1 

m 


n 

1=1 


d2, 

2'Kiz, 


'^^^lBiZa,Zc-,R)Iv{Zc)lB{Zc,Zb-,-R)Iv{Zb)=K ^ 5{Za,a{Zb)) 


(3.33) 
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Figure 8: The 46-space representations for the brane configurations in figure 4(a) and 4(b). 
The number q labeling a region indicates that there is an (N, q) 5-brane on that region. 


The appearance of Seiberg duality as an integrability condition might be unexpected 
from the point of view of lattice models. From our perspective, this is not entirely surpris¬ 
ing. The integrability is simply a consequence of the topological invariance of the relevant 
TQFT with extra dimension that has line operators, and these line operators are created 
by NS5-branes. It is well known that in the brane construction of 4d A/" = 1 gauge theories, 
Seiberg duality is realized by rearrangement of NS5-branes [63, 64]. 

4 Brane tiling models 

In the discussion of brane box models in the previous section, we found it helpful to consider 
a deformation of the brane configuration that splits the D5-branes across NS5-branes and 
creates (A^, ±1) 5-branes, as in figure 4. The bifundmantal multiplets represented by the 
horizontal and vertical arrows in the quiver originate from open strings whose ends are 
attached on the brane junctions. The original conhguration is recovered in the limit where 
the D5-brane segments are joined together and the (A^, ±1) 5-branes disappear, and these 
multiplets become massless in this limit. 

What happens if we split the D5-branes across every NS5-brane, and keep the distances 
between the segments finite? In such a situation, the horizontal and vertical multiplets are 
all massive and absent in the infrared theory. As far as renormalization group invariant 
quantities such as the supersymmetric index are concerned, these multiplets can therefore 
be ignored. On the other hand, there are now far more ways in which diagonal multiplets 
arise. To see this clearly, let us represent brane junctions by oriented dash-dot lines, which 
we call zigzag paths. The orientation of a zigzag path is such that if it is going up, then 
there is an (N, q) 5-brane on the left and {N, g -|- 1) 5-brane on the right for some q. For 
example, the two configurations in figure 4 are represented on the 46-torus as in hgure 8. 
Diagonal multiplets arise from open strings localized near intersections of rapidity lines 
or zigzag paths. When the D5-branes split, each rapidity line splits into two antiparallel 
zigzag paths. Correspondingly, each intersection point breaks apart into four intersection 
points. As the number of intersection points increases, we expect to get more diagonal 
multiplets. 

proved in [62], where Sn is the symmetric group of degree N and the delta function S is defined with respect 
to the measure f ‘tLZaj/“inizaj. The overall factor K = Ils^ii P’'l)) removed 

by rescaling the R-matrix. 
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The number of U(l) flavor symmetries also gets larger. Recall that these come from 
global gauge symmetries on the NS5-branes. In a brane box model, gauge transformations 
on the upper and lower halves of an NS5-brane give the same flavor symmetry in the 4d 
theory since they simply act in opposite manners on the horizontal and vertical multiplets 
(and the action on the diagonal multiplets is determined by the requirement that the super¬ 
potential be neutral). After the brane configuration is deformed, there are no horizontal 
or vertical multiplets anymore, so gauge transformations on the two halves can lead to 
different flavor symmetries. 

We can make these statements more precise when the deformed system contains only 
(A^, 0) and {N,±l) 5-branes. In fact, in such a situation the system is still described by 
an AA = 1 quiver gauge theory [10]. The quiver is such that an SU(A^) node is assigned to 
each (N, 0) region, and an arrow connects two nodes on regions sharing a vertex: 


i i 

I ° 

*- 1 -► ^ ^ 

0 I 1 (J I 

j j 


(An (A^, ±1) 5-brane is dual to a single D5-brane, so it supports a U(l) gauge field which 
is frozen at low energies.) To the ith zigzag path is associated a flavor symmetry U(l)i, 
and arrows intersecting with this path have nonzero U(l)i-charge T). We choose the nor¬ 
malization for the Fi such that the arrow in the above picture has Fi = +1 and Fj = — 1. 
It is known that the U(l)i account for all nonanomalous U(l) flavor symmetries present 
in the gauge theory [65]. They are not independent, however. Since each arrow is in¬ 
tersected precisely by two zigzag paths assigning opposite charges to it, the sum of the 
charges vanishes: 

i 

Hence, the number of U(l) flavor symmetries is equal to the number of zigzag paths mi¬ 
nus one. 

Once we start thinking in terms of zigzag paths, we can construct many more M = 1 
theories. We can consider various configurations of zigzag paths, not necessarily obtained 
by deforming brane box configurations. Such configurations of zigzag paths are called brane 
tilings. Provided that only {N, 0) and {N, ±1) 5-branes appear on the graph of zigzag paths, 
brane tilings yield J\f = 1 quiver gauge theories, whose quivers are determined by the same 
rule as above. Actually, brane box models can also be described as a brane tiling model; 
see figure 9 for the brane tiling configuration for the 2x3 brane box model, whose quiver 
is shown in figure 2(b). In this sense, brane tilings generalize the brane box construction. 
The reader may consult the reviews [66, 67] for more extensive discussions on brane tilings 
and their relations to other subjects such as the AdS/CFT correspondence. 

Now, let us try to connect brane tilings to an integrable lattice model. Since the 
only difference compared to the brane box case is the type of defects inserted in the 6d 
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Figure 9: The 2x3 brane box model as a brane tiling model. 
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Figure 10: The Yang-Baxter equation for zigzag paths. 


theory on the D5-branes, much of the previous argument carries over. Using this 6d theory, 
we can define a 2d TQFT whose correlation function for a configuration of zigzag paths 
is given by the supersymmetric index of the corresponding 4d AA = 1 theory. An extra 
dimension emerges if we utilize dualities and embed the 6d theory in M-theory. Zigzag 
paths are mapped to semi-inhnite M5-branes, and their x^^-coordinates provide spectral 
parameters. So we indeed have all ingredients required for the construction of an integrable 
lattice model. It is reassuring that we have found the relation (4.2), which implies that the 
supersymmetric index is invariant under an overall shift of the a:^*^-coordinates. 

The next step would be to write down an integral expression for the supersymmetric 
index and read off the R-matrix. Here we get into trouble. Clearly, we cannot do this for 
an arbitrary configuration of zigzag paths; in general, we find {N,q) regions with |g| > 1, 
and do not know what the corresponding 4d theory is.® Then we are forced to restrict our 
attention to configurations without such regions. Sadly, this restriction is not compatible 
with integrability: the Yang-Baxter equation for three zigzag paths always involves regions 
with |g| > 1, as can be seen from figure 10. 

The difficulty stems from the fact that zigzag paths are not really line operators. They 
are more properly called domain walls since they separate regions with different physical 
properties. This is a problem we did not encounter for brane box models. In that case, 
rapidity lines were genuine line operators. 

® An exception is the case with N = 2. In this case, every {N, q) 5-brane can be transformed to [N, 0) 
or (Y, 1) 5-brane by an SL(2;Z) transformation. So there are only two equivalence classes of 5-branes, 
which we can distinguish by painting black or white. This gives a checkerboard pattern on . Presumably, 
our construction applied to this case gives an integrable lattice model discovered in [1]. The Yang-Baxter 
equation for this model hrst appeared in [68], and is equivalent to an identity proved in [60]. See [69, 70] 
for more discussions and further developments in this direction. 
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Our consideration at the beginning of this section suggests a resolution to this problem. 
We have seen that when a brane box model is deformed to a brane tiling model, a rapidity 
line splits into two antiparallel zigzag paths. Conversely, we can take two neighboring 
antiparallel zigzag paths, and regard them as one thick line (or “ribbon”) operator. In this 
way we can make two kinds of line operators, namely 


(ri.ra) 



(4.3) 


and another one for which the orientations of the zigzag paths on the right-hand side are 
flipped. Below we will only use the first one. 

If we use this line operator and give the “background charge” g = — 1 to the TQFT, 
then regions with |g| > 1 never arise: 



(4.4) 


Thus, a brane tiling constructed with line operators of this type is always described by a 
quiver gauge theory, and we can readily write down an integral expression for its supersym¬ 
metric index. From this expression we can read off the R-matrix of the associated vertex 
model. Just as in the case for brane box models, we expect that the Yang-Baxter equation 
for this R-matrix can be understood as an equality between the supersymmetric indices of 
two dual gauge theories. Let us see if this is the case. 

From the quiver rule (4.1), we hnd that the lattice model in question is a vertex model^® 
with the R-matrix given by 


32 

*2 



(si, 

S2) 



(4.5) 


Figure 11 shows the case when the model is placed on a 2 x 3 lattice. The R-charges of the 
bifundamental multiplets, specified by the numbers accompanying the arrows, are chosen 
to be consistent with the normalization of the flavor charges and add up to 2 around loops. 
With this choice, the R-charges satisfy the anomaly cancellation condition. 

What we have just constructed is the integrable lattice model studied in [2, 8]. (Its 
generalization to the lens space index was considered in [9].) The Yang-Baxter equation 

^°The same model can also be described as an IRF model. However, the IRF description uses line 
operators of both types, and the Yang-Baxter equation involves twisting of one of the three ribbons. See [9] 
for a discussion from the IRF viewpoint. 
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(a) 



(b) 


Figure 11: (a) The 2x3 lattice model constructed from the line operator (4.3). (b) The 
quiver of the corresponding gauge theory. 


for the model is 



For simplicity we have omitted the R-charges in the quivers. As expected, this equation 
follows from Seiberg duality: applying the basic duality transformation (3.31) four times, 
we can turn the left-hand side into the right-hand side. This process is illustrated in 
figure 12. 

Incidentally, a single Seiberg duality transformation is induced by the Yang-Baxter 
move involving two zigzag paths and one “twisted” ribbon operator: 



(4.7) 


This relation between the Yang-Baxter move and Seiberg duality was pointed out by 
Hanany and Vegh [71]. 



5 3d A/” = 2 and 2d M = (2,2) quiver gauge theories 

The results from the previous sections imply, via T-duality, that there exist similar con¬ 
nections between lower-dimensional quiver gauge theories and integrable lattice models. 
In this section we discuss integrable lattice models associated with 3d AA = 2 and 2d 
N = (2, 2) theories. We will focus on those coming from brane box configurations, though 
our analysis can be adapted straightforwardly to the brane tiling case. 
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Figure 12: A sequence of duality transformations that demonstrates the Yang-Baxter 
equation (4.6). The letter on an arrow between two quivers indicates the node at which 
the transformation is applied. Note that circle nodes can be relabeled freely since their 
fugacities are integrated over. 


5.1 3d brane box models 

Let us consider the following brane setup in Type IIA superstring theory, which is related 
to the 4d brane box model (3.1) by T-duality along the 2-direction: 
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This is a brane box configuration that realizes a 3d AA = 2 quiver gauge theory. Being re¬ 
lated by T-duality, the theory is essentially the dimensional reduction of the corresponding 
4d AA = 1 theory, and is described by the same quiver (with all Chern-Simons levels van¬ 
ishing). We expect that the supersymmetric index of this theory is given by the partition 
function of an integrable lattice model. 

As in the 4d case, the theory can be placed on 5^ x with all four supercharges 
unbroken [72, 73]. We are interested in the supersymmetric index 

X(x; {r,}) = TV 52 ((- 1 )^x'^+ 2 ^ 3 ^ ^ ( 5 . 2 ) 

where J 3 is the Cartan generator of the isometry group SU(2) of 5^. It receives contribu¬ 
tions only from states with energy E = R J 3 . As usual, we have omitted the fugacities 
Ua for the flavor symmetries U(l)a, whose effects can be taken into account by keeping the 
dependence of the R-charge R on the spectral parameters Tq,. 
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The 3d supersymmetric index is computed by an integral formula similar to the one in 
the 4d case. The novelty is that there can be nontrivial magnetic fluxes through S‘^. For 
each gauge group SU(A^)a, the flux is specified by a set of integers rua = , ma,N) 

subject to the condition rua^i = 0. Therefore, the formula involves a sum over fluxes as 
well as an integration over the fugacities for the gauge groups: 


1 = 



N-l 


n 


27riZaj 




rUa) 


n ,ma), {zb,m!,);Rab). 
©-►® 


(5.3) 


Here | Wma I is the order of the Weyl group of the subgroup of G that is left unbroken by 
rua, and the vector and bifundamental multiplet factors are given by [72, 74, 75] 


lY{Za,ma) = (5 4 ) 


ZB(XZa,ma),{zb,mb)]Rab) =W{x^ 






Ix‘^)oo 

The bifundamental factor Ib satisfies the identity 

ZB{{Za,ma),{zb,mb)-,Rah)l'B{{zb,mb),{Za,ma)]Rha) = 1 


(5.5) 


(5.6) 


for Rab T — 2. 

We can connect the supersymmetric index to a lattice model by following the same 
procedure as before. We define an IRF model on the lattice drawn on the 46-torus by 
the NS5-branes. The state variables are {za,ma), and the Bolzmann weight is given by 
the same picture (3.27). By construction, the supersymmetric index of a 3d brane box 
model coincides with the partition function of this IRF model. The Yang-Baxter equation 
again boils down to the equality (3.31) between the supersymmetric indices of two gauge 
theories. This equality, if true, establishes the integrability of the IRF model. 

The equality indeed follows from the SU(W) version [76, 77] of the duality for 3d 
M = 2 gauge theories proposed by Aharony [19]. There is an interesting twist, though: the 
theories described by the quivers that appear in the equality are not dual to each other. 

The correct duality takes the theory on the left-hand side to a theory whose matter 
content consists of the chiral multiplets described by the quiver on the right-hand side, 
plus two extra singlets 17+, V-. More generally, Aharony duality relates two 3d A/" = 2 
gauge theories with the following matter contents: 
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The gauge group G is SU(A^c) for the theory with Q, Q, and SU(A^/ — Nc) for its dual. 
Although a proof has not been given yet, there is strong evidence for the equality of the 
supersymmetric indices of these two theories in terms of series expansion in x [77-79].^^ 

The above table contains one more global symmetry compared to its 4d counter¬ 
part (3.32), namely the axial symmetry U(l)y 4 which is anomalous in four dimensions. 
In the dehnition of the index (5.2), we did not include fugacities for axial symmetries. 
Thanks to this omission, V± always have R = 1 for whatever shift of R by the relevant 
U(l) flavor symmetries, and their contributions to the index cancel since we can turn on 
a superpotential W = mV+V-. If we included these fugacities, the cancellation would no 
longer occur. In that case, the presence of V± would be crucial for the indices of the two 
theories to match. 

Why does the logic fail if the fugacities for axial symmetries are turned on? The 
point is that when we apply T-duality to get a 3d brane box model from a 4d one, the 
2-direction is compactified to a circle, however small its radius may be. In other words, we 
are really compactifying the 4d theory on with finite radius, not dimensional reducing 
it. As emphasized in [76], the effective 3d theory obtained by compactification has a 
superpotential that contains monopole operators. This superpotential breaks the axial 
symmetries, thereby reproducing the effect of anomalies in the 4d theory. So it is not that 
the logic fails. Rather, we have no fugacities to turn on to begin with. Upon inclusion of 
such superpotentials on both sides, the extra multiplets disappear from the duality. 

The axial symmetries are recovered if we decompactify the 2-direction. Although it is 
perfectly hne to do so, the resulting system is not related to a 4d brane box model by T- 
duality and its connection to a lattice model requires a more careful analysis. It is not clear 
whether these additional symmetries admit a natural interpretation in the lattice model. 

5.2 2d brane box models 

Applying further T-duality along the 1-direction, we obtain the following brane configura¬ 
tion realizing a 2d A? = (2, 2) quiver gauge theory: 
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This theory is described by the same quiver diagram as in the 3d and 4d brane box cases. 
However, each node now represents a U(AI) gauge group, not SU(A^). The reason is that 
unlike the case with D4- or D5-branes, motion of D3-brane segments along NS5-branes 
creates disturbances on the NS5-branes that fall off sufficiently rapidly, and costs only a 
finite amount of energy. (In the 3d brane box case, on the other hand, the required energy 
becomes arbitrarily large as we increase the radius of the 2-circle.) 

^^For Nc = 2 and Nf = 3, 4, the equality was recently proved in [80] in the absence of the fugacity for 
U(1)a. 
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The supersymmetric index of a 2d A/" = (2, 2) theory is known as the elliptic genus. It 
is given by 

Aq, y; {’’«}) = Tr^i , (5.9) 

where is the left-moving Hamiltonian, and R is defined as twice the left-moving R- 
charge so that a superpotential preserving the R-symmetry has R = 2. The trace receives 
contributions only from states with the right-moving Hamiltonian Hr = 0. The elliptic 
genus can be computed as the partition function with twisted boundary conditions on 
a torus whose complex parameter r is given hy q = As such, it enjoys a nice 

modular property. 

An integral formula for the elliptic genus of an AA = (2, 2) gauge theory was derived 
in [81-83]. For the theory we are considering, the formula reads^^ 

r 1 ^ d 

^ = f n M n ^^v(na) H lBiUa,Ub-,Rab), (5.10) 

@ 1=1 @-*-® 

(5.11) 

(5.12) 

Here z is related to the fugacity y hy y = and 


with the vector and bifundamental multplet factors given by 
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PljUg^i - Uqjlr) 
t 0liUa,i - UaJ - z\t)'' 


Rab) — 




^li'^bj '^a,i “ 1 “ RabZ/‘^\T^ 


r?(r)=gV24n(l-g-), 

n=l 


0i{z\t) 


OO 

^ ^_^'^n^2TTiz{n+l/2)^7riT{n+l/2)^ 

n=—oo 


(5.13) 

(5.14) 


are the Dedekind eta function and a Jacobi theta function. Note that 6i{—z\t) = —0i{z\t). 
As a consequence, Xb obeys the identity (3.28). 

The above formula takes the same form as the 4d index formula (3.24) under the 
identification Za^i = The subtlety lies in the integration contour, as was elucidated 

in [82, 83] with a careful path integral analysis. Since 0i{z\t) has a pole at z = 0, the 
integrand has poles at various places in the u-space, and the contour must pick up the 
residues at some but not all of these poles. There are many choices that lead to the correct 
answer. For us, a simple choice is such that for each node a, every Ua,i encircles the pole 
located at Ua,i = Ubj — Rbaz/2 for some (6, j), which comes from the factor X^iub, Ua \ Rba) 
associated to the incoming arrow 

©-- ® ■ (5-15) 

^^For concreteness we will follow the treatment in [82, 83] where the elliptic genus is defined in the (R, R) 
sector, though the considerations given below equally apply to the other sectors. See [81] for the formula 
for the (NS, NS) sector. 
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Moreover, Ua,i and Uaj for i ^ j must encircle different poles. We sum over the contribu¬ 
tions from all such contours. Alternatively, we may choose the contours to encircle poles 
coming from outgoing arrows. The two choices give the same result, up to an overall sign. 
We refer the reader to [82, 83] for more details on the integration contour in general, and 
section 4.6 of [83] for details specific to theories closely related to ours. 

As in the higher-dimensional cases, the supersymmetric index of a 2d brane box model 
coincides with the partition function of an IRF model. The Boltzmann weight is given 
by the formula (3.27), and the Yang-Baxter equation reduces to the equality (3.31). This 
time, let us demonstrate this equality explicitly. 

The elliptic genus of the theory described by the quiver on the left-hand side is com¬ 
puted by the contour integral 

^hllS ~ ^ Ug] Rag)^Biudj Ug'i Rdg) 

i=i 

X lB{ug,Uc;Rgc)lB{ug,Uf;Rgf). (5.16) 

We choose to evaluate it with contours that encircle N poles from the factors associated 
with the incoming arrows. For instance, suppose that we are evaluating the integral for 
contours such that the unordered set of these poles is {uaj — Ragzf2}, that is, each Ug^i- 
integral picks up the residue at Ug^i = Ua,j — Ragzj2 for some j. The N\ different ways to 
assign the poles to the Ug^i make equal contributions, so we only need to consider the case 
when Ug^i = Ua^i — Ragz/2 for all i. Using the fact that the residue of 1/9i{u\t) at u = 0 is 
l/27rr/(r)^, we find that the sum of the contributions from these contours is 

RBiUai Uc] Rag T Rgc)RB{Ua, Uf] Rag T Rgf')'I'B{!^diUa\Rdg Rag)' (^'l"^) 

On the other hand, the elliptic genus of the right-hand side is given by 

^IrHS ~ 2^ • ^'vi^g)^B {Ug, Ud'i Rgd)RB{Ug, Ua'i Rga) 

X Xg (Uci Ug , Rcg)R'Q ) Ug, Rf g) 

^ Xg {v^a-i ^c; Rac^R-B Raf)RB '^ci RdcjR-B Rdf) • (5.18) 

For the dual theory, we choose the poles from the outgoing arrows. The contours picking 
up the residues at the poles {udf + Rgd/2} contribute 

'R'Bifdd'! ^aj Rga Rgd)RBi^ci Reg T Rgd)RBi^f ^ ^dj Rfg T Rgd) 

X Xg (Mqi , Uc, 7?ac)Tg (Ua, Uf, Ra/)Tg {udj Uc, Rdc)RB (^dj Uf, Rdf)' (5.19) 

After canceling out some factors using the identity (3.28), we see that this is equal to the 
contribution (5.17) found above. 

In a similar fashion, one can show that for each set of poles {%,i} C {uaf — Ragzj2} U 
{udf — Rdgzf2} in the theory on the left-hand side, there is a set of poles {ug^i} C {uaf + 
Rgazj2} U {ud^i + Rgdzj2} in the theory on the right-hand side such that the contours 
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encircling these poles give the same contributions to the elliptic genera of the respective 
theories. (In the case that the R-charges vanish, {%,i} is the complement of {ug,i} in 
{ua^i} U {ud^i}.) This is a one-to-one correspondence between the choices of poles in the 
two theories. Therefore, Xlhs = 2iRHS- 

The equality just proved is basically a consequence of a variant of Hori-Tong dual¬ 
ity [20] proposed in [84]. However, it is important that the fugacities Ua, Uc, u^, uj are 
unconstrained here. By contrast, when regarded as the elliptic genera of dual theories, the 
two sides of the equality must be evaluated with the constraint Ylii^a,i+Uc,i+Udd+Uf^i) = 0 
since the diagonal U(l) subgroup of the flavor group U(A^)a x \J{N)c x \J{N)d x U(A^)/ is 
gauged. So strictly speaking, the duality does not imply the Yang-Baxter equation. Also, 
note that the fugacity for the axial symmetry U(l)^ is set to zero as in the 3d case. In 
general, this can be turned on in the elliptic genus. 

Finally, we remark that the duality does not work if we replace the gauge group with 
S\]{N). This is unlike the higher-dimensional cases, where one can gauge U(1 )b to obtain 
the duality for V{N) theories from that for SU(A^) theories. Had we started with the 
dimensional reduction of the 3d brane box model and simply discarded the U(l)^ fugacity, 
we would have reached a wrong conclusion. 

6 2d Af = (0, 2) quiver gauge theories 

Looking at how the duality relation (3.31) was verified for 2d brane box models, we notice 
that the numerator and denominator of the bifundamental factor Xb played rather different 
roles. More specifically, while the denominator determined the pole structure and hence 
the possible choices for integration contours, the numerator only provided cancellations 
of some factors. The separation in their roles hints at the existence of integrability in 
less supersymmetric situations where an M = (2, 2) chiral multiplet decomposes into two 
multiplets that correspond to the numerator and denominator. 

M = (0, 2) supersymmetry provides precisely the required decomposition: an M = 
(2, 2) chiral multiplet consists of an AA = (0, 2) chiral multiplet and Fermi multiplet, with 
the former corresponding to the denominator and the latter to the numerator. In this 
section we introduce three classes oi N = (0, 2) quiver gauge theories whose elliptic genera 
are captured by integrable lattice models. The first two classes lead to IRF models much 
like N = (2, 2) brane box and brane tiling models discussed in section 5.2. The third class 
gives rise to an IRC model on 3d lattices. 

6.1 Af = (0, 2) theories related to brane box configurations 

Let us take a quiver from some brane box model, and replace the vertical arrows with dotted 
arrows. (Of course, one may as well choose to replace the horizontal arrows instead.) For 
example, starting from the quiver for the 2x3 brane box model shown in figure 2(b), 
we obtain the diagram in figure 13. We interpret the resulting diagram as a quiver of an 
N = (0, 2) gauge theory, letting dotted arrows represent bifundamental Fermi multiplets. 
As usual, circle nodes and solid arrows represent vector multiplets and bifundamental chiral 
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Figure 13: The M = (0, 2) theory obtained from the 2x3 brane box model. 


multiplets. By this procedure we can associate an A/" = (0, 2) quiver gauge theory to every 
brane box model. The reader is referred to [21, 85] for basics of Af = (0, 2) gauge theories. 

Although the left-moving R-symmetry is no longer present, the theory still has the 
flavor symmetries U(l)ci, at least classically. So we define the “R-charge” in the present 
case by 

R = '^raFa. (6.1) 

a 

Our claim is that the elliptic genus (5.9) of this theory is equal to the partition function of 
an integrable IRF model. 

Before we give a proof of this claim, cautionary remarks are in order. As it is, the 
theory defined above suffers from mixed anomalies for the diagonal U(l) factors of the 
U(A^) gauge groups. Similarly, the flavor symmetry generated by R is anomalous. These 
anomalies must be canceled by introduction of extra multiplets. For the moment let us 
ignore these issues; we will address them later. 

The relevant formula for elliptic genera was derived in [82, 83]. The only difference 
compared to the Af = (2, 2) case is that it involves three kinds of factors, corresponding to 
the three kinds of multiplets: 

^ ^ / n n n ^'B.c{Ua,Uh]Rah) F.BF{Ua,Uh]Rah), (6.2) 

@ i=i @-^® (a) I-® 


where the vector, bifundamental chiral and bifundamental Fermi multiple! factors are^^ 


Tv(ua) 


TBc(^^a, Uh] Rab) 


F-BYiUa-, Ub'i Rab) 


(^) n- 

n 


i0liua,i - Uajlr) 


iAj 

iriir) 


r]{T) 




n 


'^a,i T RabZ/‘^I'r) 
iOxi^UbJ T Rab^/‘^l’^) 


r]{T) 


(6.3) 

(6.4) 

(6.5) 


Up to an immaterial overall sign, I^dua, Ub', Rab)FBF{ua, Ub] Rab — 2) is equal to the Af = 
(2, 2) bifundamental chiral multiple! factor (5.12). Likewise, Iy{ua)FBciua, Ua] —2) is equal 
to the Af = (2,2) vector multiple! factor (5.11), as is consistent with the decomposition 

^^See [21, 86] for the formulas for the (NS, NS) sector. 
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of an AA = (2, 2) vector multiplet into an AA = (0, 2) vector multiplet and adjoint chiral 
mnltiplet. We have 

1-Bc{Ua,U}y] Rab)1-BF{Ub-,Ua] Rba) = 1 (6.6) 

for Rab + Rba = 0, or more graphically, 


HCS = E B ■ 

The cancellation of arrows corresponds to giving masses to a pair of chiral multiplet and 
Fermi multiplet. 

For the theory under consideration, the elliptic genus can be computed as the partition 
function of an IRF model with Boltzmann weight 


r 



s 



( 6 . 8 ) 


Plugging this into the Yang-Baxter equation, we hnd that the integrability of the model is 
equivalent to the equality 



Let us demonstrate it. 

The left-hand side is given by the contour integral 


Xlhs 



27ri 


^v(^g)^Bc(^(i) '^gi Rdg)^TiF{Uai Ug, Rag) 

X '^BciUg, Uc] Rgc)'^Bc{Ug,Uf; Rgf)- (6.10) 


If we choose the contours to pick up residues from the incoming arrow, then the set of 
relevant poles is {ud,* ~ Rdgzl2}. Thus we have 

2 lHS = ^BF iUa,Ud-, Rag - Rdg)^Bc{Ud, Uc, Rdg + Rgc)^Bc{Ud, Uf, Rdg + Rgf) ■ (6.11) 

The right-hand side of the equality (6.9) is 


2^rhs 


1 

m 


N 


dUn 


i=i 


X Tbc(^C5 Ug, Rag)RBc{fdif t Ug, Rfg) 

X lBF{Ua, Uc] Rac)lBF{Ua, Uf] Raf)RBc{Ud, Uc] Rdc)'^BciUd, Uf] Rdf)- (6.12) 
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Figure 14: 


The R-charges of the extra multiplets for anomaly cancellation. 


The “complement” set of poles is {ua,i + Rgazl2}, and they give 


^RHS — TBp(?rQ,, ltd, ^gd R(jia)TBc(^C) Ua, Reg T Rga)^Bc{^f, ^fg T Rga) 

X XBF(Wa, Uc, Rac)TBF('Ua, Uf] Raf)RBciUd, Uc', Rdc)^BciUd, Uf, Rdf)- (6.13) 


Canceling some factors using the identity (6.6), we find Xlhs = Xrhs- 

In fact, the quivers on the two sides of the equality (6.9) describe two theories related 
by the (0, 2) triality transformation [21]. The third theory in the triality is the free theory 
with quiver 



(6.14) 


This is manifest in the expression (6.11) of the elliptic genus. 

Finally, we come back to the issues of anomalies. To cancel the mixed U(l) gauge 
anomalies and R-anomaly, we add to each solid arrow 


(a)-(6-15) 

a Fermi multiplet with R = —NRab—ein the representation deta(8)det^^ of U(A^)aXU(A^)b, 
and to each dotted arrow 

0 (6.16) 

a chiral multiplet with R = —NR^b + e in deta <8> det^^. Here e is a nonzero parameter. 
Furthermore, to each node we introduce a singlet chiral multiplet with R = e, though this 
is not required for anomaly cancellation. We can visualize these extra multiplets as wavy 
arrows, as in figure 14. 

The presence of the extra multiplets modifies the factors associated to vector and 
bifundamental multiplets to 




^BF Ub, Rab) 


r/(r) 

ir]{T) 


^l(Ei ^a,i - Ej Ub,i + {-NRab + e)z/2\T) 


1-BF{Ua, Ub] Rab)- 


(6.17) 

(6.18) 
(6.19) 
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Figure 15: The M = (0, 2) theory obtained from the 2x3 brane tiling model. 


Apart from this modification, the computations of the elliptic genera for the equality (6.9) 
remain unchanged. In particular, we can use the same integration contours, since the extra 
multiplets introduce neither solid incoming arrows to the left-hand side nor solid outgoing 
arrows to the right-hand side. 


6.2 Af = (0, 2) theories related to brane tiling configurations 

Quivers from brane tiling models can also be modified to produce Af = (0, 2) theories that 
exhibit integrability. As in the brane box case discussed above, we obtain these theories 
by replacing one of the arrows in the R-matrix with a dotted arrow: 



(S1,S2) 


( 6 . 20 ) 


Extra multiplets for anomaly cancellations are implicit in this picture. Figure 15 shows the 
quiver for the Af = (0, 2) theory obtained from the 2x3 brane tiling model, whose quiver 
is given in figure 11(b). 

The Yang-Baxter equation for this R-matrix reduces to the equality 




( 6 . 21 ) 


Application of the (0, 2) triality transformation (6.9) four times (more precisely, one triality 
transformation followed by three others in the opposite direction) turns the left-hand side 
into the right-hand side. This is illustrated in figure 16. 


6.3 Brane cube models 

Lastly, we briefly discuss the class oi Af = (0, 2) quiver gauge theories known as brane 
cube models [22]. These theories are constructed from brane configurations analogous to 
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Figure 16: A sequence of (0, 2) triality transformations that relates the two sides of the 
Yang-Baxter equation (6.21). An arrow between two quivers represents a (0,2) triality 
transformation in the direction from the right- to left-hand side of the equality (6.9). 


brane boxes: 
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Note that brane cube models include 3d brane box models (5.1) as special cases. A brane 
box model is characterized by the 2d lattice drawn on the 46-torus by the NS5-branes. For 
a brane cube model, the NS5-branes make a 3d lattice in the 246-space which we take to 
be a 3-torus T^. 

A brane cube conhguration can be thought of as describing intersecting codimension-1 
defects in the 5d theory on the D4-branes. We can topologically twist the theory along the 
246-space, replacing the holonomy group 80(3)246 with the diagonal subgroup of 80(3)246 x 
80 ( 3 ) 357 . The twisted theory has four supercharges that are scalar in the 246-space, hence 
unbroken when this space is a general 3-manifold S. These supercharges generate J\f = 
(2, 2) supersymmetry in the 01-space. Half of them, generating an N = (0, 2) subalgebra, 
are preserved by codimension-1 defects supported on arbitrary surfaces inside S. Therefore, 
the twisted theory gives a morphism that sends a configuration of surfaces in S to the 
effective M = (0, 2) theory. Taking the elliptic genus, we get a 3d TQFT on S equipped 
with surface operators. An extra dimension emerges upon lifting the system to M-theory, 
and it follows that the elliptic genus of a brane cube model is equal to the partition function 
of a 3d integrable lattice model. 

Like an IRF model is associated with brane box models, the lattice model associated 
with brane cube models is an IRC model. According to the quiver rule found in [22], one 
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way to describe this model is to set its Boltzmann weight 



(6.23) 


This Boltzmann weight, or another one that leads to the same IRC model on cubic lat¬ 
tices, should solve the tetrahedron equation (2.25). It would be interesting to verify this 
statement. 
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